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Allocation of limiting resources, such as nutrients, is an important adaptation strategy for plants. Plants 
may allocate different nutrients within a specific organ or the same nutrient among different organs. In this 
study, we investigated the allocation strategies of nitrogen (N) and phosphorus (P) in leaves, stems and roots 
of 126 shrub species from 172 shrubland communities in Northern China using scaling analyses. Results 
showed that N and P have different scaling relationships among plant organs. The scaling relationships of N 
concentration across different plant organs tended to be allometric between leaves and non-leaf organs, and 
isometric between non-leaf organs. Whilst the scaling relationships of P concentration tended to be 
allometric between roots and non-root organs, and isometric between non-root organs. In arid 
environments, plant tend to have higher nutrient concentration in leaves at given root or stem nutrient 
concentration. Evolutionary history affected the scaling relationships of N concentration slightly, but not 
affected those of P concentration. Despite fairly consistent nutrients allocation strategies existed in 
independently evolving lineages, evolutionary history and environments still led to variations on these 
strategies. 

Allocation of limiting resources is an important strategy for plants to adapt to changes in their envir- 
onment 1 , which has been widely observed in the allocation of biomass 2 4 , functional traits 5 and morpho- 
logical plasticity 6,7 . Nutrients availability plays an important role in ecosystem function and development, 
since some essential nutrient elements, such as nitrogen (N) and phosphorus (P), are considered limiting 
resources in nature 8,9 . Allocation of nutrients involves in two main ways: strategic allocation of different nutrients 
within a specific organ and allocation of the same nutrient among organs. It has been widely observed that how 
different nutrients are coupled within a specific organ, such as leaves 10 " 12 , stems 13 , fine roots 14,15 and twigs 14 . 
However, relatively little is known about how the same nutrient is allocated among organs. Plants need to allocate 
nutrients they absorb into different organs to meet the requirements of multiple functions such as growth, 
reproduction, nutrient storage and defense 1 . Nutrient allocation also reflects a plant's relative investments of 
essential molecules to different organs, such as N-rich proteins and P-rich RNAs 16 . 

Previous studies have demonstrated that evolutionary history, environmental stresses and plant functional 
groups are possible factors influencing the relationship of nutrients among different organs 16 " 19 . From an evolu- 
tionary perspective, closely related species may share a similar evolutionary history and contain life history 
strategy 5,20,21 . Exploring how plants with different evolutionary history allocate N and P among organs can 
improve our understanding of correlated evolution 16 . In a meta-analysis, Kerkhoff et al. 16 found that the scaling 
relationships of N and P concentrations between plant organs changed between metabolic and structural organs 
and proposed a general relationship describing N and P allocation among organs across diverse plant lineages. 
From a life-history perspective, exploring how plants allocate N and P among organs can improve our under- 
standing of nutrient utilization strategies under climate change scenarios 22 . Plants need to balance the allocation 
of limiting resource to maximize their growth and may change allocation strategies as a response to changes in 
their environments such as climate and soil nutrient availability 1,17,18,22,24 . For example, Sardans & Penuelas 22 
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found a higher allocation of nutrients to leaves than to woody bio- 
mass in more humid climatic conditions. Plants may contain higher 
N concentrations in leaves to better adapt to arid environments 25-29 
through exploiting greater light availability 25 while reducing stomatal 
conductance and increasing water use efficiency 28-2 '. Plant functional 
groups have also been shown to influence allocation of nutrients 
through different photosynthetic pathways, fraction of woody tissue 
and N-fixation ability 16-19 . Further studies are needed for a better 
understanding of the combined impacts of climate, soil and evolu- 
tionary history on the allocation of nutrient among plant organs. 

Nutrients such as N and P frequently limit plant growth and play 
important roles in plant functioning 14,3031 . As a major element in 
proteins, N in plant tissues is related to various functions, such as 
photosynthesis in leaves, tissue respiration, conversion and storage 
of photosynthates and chemical defense in all plant organs 32,33 . As a 
major element of RNAs, P is mainly involved in the production of 
proteins 32,34 . Since different plant organs perform different functions 
in plant growth, the requirements of N and P will differ. For example, 
leaf N and P concentrations were found to have different scaling 
relationships with leaf photosynthetic capacity 23 . These differences 
may lead to differences in the scaling relationships of N and P con- 
centrations among plant organs. Additionally, N and P concentra- 
tions in plants are widely reported to be more similar among closely 
related species 35 . Finally, plants are able to coordinate various organs 
with respect to absorption and allocation of limiting resources and 
adaptation to environmental constraints 18 ; for example, plants tend 
to increase leaf N concentration to maintain growth under dry con- 
ditions 2529 . Given the aforementioned ability of plants to adapt to 
different allocation strategies, we made the three following hypo- 
theses. First, N and P have different scaling relationships across 
organs because of their different physiological functions. Second, 
the scaling relationships of N and P across organs depend on the 
phylogenetic relatedness of the plants; more similar nutrient alloca- 
tion strategies exist in more closely related species. Third, we hypo- 
thesize that the scaling relationship of nutrients among plant organs 
changes with environmental factors, such as soil nutrients and pre- 
cipitation; particularly, we expect higher leaf N concentration at 
given root or stem N concentration in more arid regions. 

We tested the aforementioned hypotheses using an extensive 
investigation of N and P storage in different organs of shrubs across 
northern China. Compared to trees, shrubs are generally smaller and 
relatively more uniform in size among species because they do not 
have large trunks or roots, which weaken the "dilution effects" of N 
and P allocation to structural components of woody plants 16 . 
Shrubland represent a major vegetation type in Northern China with 
community types varying from east to west China along an aridity 
gradient 36 . In this study, we explored the relationships of N and P 
among leaves, stems and roots of shrubs in Northern China using 
scaling approach, which has been applied in exploring generalities of 
patterns and trade-offs in many plant traits 10,16 . Nutrient allocation 
among plant organs involves two facets: concentration and biomass 
storage 16 . We compared the scaling relationship for both concentra- 
tions and storages of N or P among leaves, stems and roots. The 
comparison has been conducted at following scales: 1) individual vs. 
species vs. phylogenetic levels; 2) species from the family Fabaceae 
vs. those from other families; 3) samples from different soil nutrient 
concentration levels and 4) samples from different climatic regions. 

Results 

Scaling of N and P among organs. Both concentrations and storages 
of N and P showed significant correlation among all organs at the 
individual, species, and phylogeny levels (Table 1). 

At the individual level, the slopes of N concentration in all organ 
pairs were significantly different from 1; the slopes for root vs. leaf 
and stem vs. leaf were significantly larger than that of root vs. stem. 
At the species level, the slopes were all significantly different from 1; 
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Figure 1 | Scatterplots showing the RMA regressions of N concentration (a-c), P concentration (d-f), N storage (g-i) and P storage (j-1) among 
organs for legume (black solid dots) and non-legume (gray open cycles) shrubs. Separate lines indicate the slopes of non-legumes (blue) and legumes 
(red) are significantly different (likelihood ratio tests, P < 0.05), whereas a single black line indicates otherwise. Lines with the slopes equal to 1 are shown 
with dotted lines. 



the slope was significantly larger for root vs. leaf and stem vs. leaf 
than for root vs. stem. At the phylogeny level, the PIC analyses also 
revealed a similar pattern to individual and species level (Table 1). 
The slope of the root vs. stem was not significantly different from 1 
when PIC is applied. There were no significant differences between 
the slopes at individual level and species level (P > 0.05). Slopes at the 
species level were not significantly different (P > 0.05) from the 
phylogeny level except the root vs. stem (P = 0.012). 

Scaling relationship of P concentration among organs was differ- 
ent from that of N (Table 1). At the individual level, the slopes of P 
concentration in all organ pairs were significantly different from 1; 
the slope for root vs. leaf or root vs. stem was more deviated than for 
stem vs. leaf. At the species level, the slope was significantly larger for 
root vs. leaf and root vs. stem than for stem vs. leaf (Table 1). The PIC 
analysis depicted that, at the phylogeny level, the slope was signifi- 
cantly larger for root vs. leaf and root vs. stem than for stem vs. leaf. 



There were no significant differences between the individual, species 
and phylogeny level (P > 0.05). 

The RMA slopes for N and P storages of the same organ pairs were 
similar at all three levels. For N storage at individual and species 
levels, the slopes of root vs. leaf, stem vs. leaf and root vs. stem were 
significantly larger than 1. The PIC analysis depicted that, at the 
phylogeny level, the slopes of root vs. leaf and root vs. stem were 
significantly larger than 1, while that of stem vs. leaf was equivalent to 
1 . For P storage at all levels, the slopes of root vs. leaf and root vs. stem 
were significantly larger than 1, while the slopes of stem vs. leaf were 
equivalent to 1. 

Influence of functional type on scaling of N and P among organs. 

All RMA regression slopes were not significantly different between 
legumes and non-legumes (P > 0.05), except the root N concen- 
tration vs. stem N concentration (P = 0.04) (Figure 1). 



SCIENTIFIC REPORTS | 4:5448 | DOI: 10.1038/srep05448 



3 



All RMA regression slopes of N and P storages among different 
plant organs were not significantly different between legumes and 
non-legumes (P > 0.05). 

Influence of soil nutrient on scaling of N and P among organs. The 

slopes of N concentration among organs were not significantly 
different across different soil total N concentration (STN) levels, 
except that at the highest STN level. The slope of root vs. stem was 
significantly different than the slopes at the other two STN levels 
(high STN vs. low STN: P = 0.009 and high STN vs. medium STN: P 
= 0.005) (Figure 2a). The slopes of N storage among organs were not 
significantly different across different STN levels (P > 0.05) 
(Figure 2d). 

At different soil total P concentration (STP) levels, the slopes of P 
concentration among organs were not significantly different across 
different STP levels, except a slightly larger slope of root vs. stem at 
the highest STP (high STP vs. medium STP: P = 0.017) (Figure 2b). 
The slopes of P storage among organs were not significantly different 
across different STP levels (P > 0.05) (Figure 2e). 

Influence of aridity on scaling of N and P among organs. The 

slopes of N concentration in root vs. leaf were significantly higher 
in semi-arid/semi-humid (Semi) and humid (Humid) regions than 
in arid regions (Arid) (Arid vs. Semi: P = 0.020 and Arid vs. Humid: 
P = 0.014). Similarly, the slopes of stem vs. leaf N concentration were 
significantly higher in Semi and Humid than in Arid (Arid vs. Semi: 
P = 0.031 and Arid vs. Humid: P = 0.049). The slope of root vs. 
stem was not significantly different among three climatic regions 
(Figure 2c). 

The slope of P concentration in root vs. stem was not significantly 
different between Arid and Semi, and significantly higher in Humid 
than in Semi (Arid vs. Humid: P = 0.048 and Semi vs. Humid: P < 
0.001). The slopes of stem vs. leaf and root vs. stem were not signifi- 
cantly different among three climatic regions (Figure 2c). 

The slopes of N and P storages among different organs were quite 
similar at different climatic regions (Figure 2f). The slopes of N and P 
storage in root vs. leaf were not significantly different between Arid 
and Semi, but were significantly lower in Humid than in Semi (P = 
0.001 and 0.006 for N and P, respectively). The slopes of P storage in 
stem vs. leaf were significantly higher in Arid than in Semi (P = 
0.028), and were not significantly different between Semi and 
Humid. Slopes of N storage were not significantly different among 
three climatic regions. The slopes in root vs. stem were significantly 
higher in Semi than in Arid and Humid (Arid vs. Semi: P = 0.001 for 
both N and P; Semi vs. Humid: P < 0.001 for both N and P). 

Discussion 

Using concentrations and storages of N and P in leaves, stems and 
roots of 126 shrub species from 172 shrubland sites, we studied the 
allocation strategies of N and P among organs through a scaling 
approach. Both concentrations and storages of N and P between 
all pairs of plant organs are significantly correlated at individual, 
species and phylogeny levels (Table 1). The coordinate variations 
of nutrients across different organs demonstrate that plant organs 
are not independent. The allocation of nutrients among organs might 
be another important part of life history strategy for plants. For 
example, plants with high leaf nutrient concentration would also 
have high nutrient concentration in stem and root to increase nutri- 
ent uptaking and phloem loading to meet the requirement of higher 
photosynthesis rate and photosynthate export. At the community 
level, this coordination would also help to predict changes in nutrient 
storage in plants under the climate change scenarios, since increased 
dominance by species with higher nutrient concentration and stor- 
age in one organ will lead to predictable increases in nutrient con- 
centration and storage in other organs 16 . 



Concentrations of nutrients in lower organs increase with those in 
the neighboring upper organs in a power law form (e.g., roots vs. 
stems, stems vs. leaves, see Figures 1 and 2, Table 1). This result is 
consistent with Brouwer's hypothesis that plant organs are compet- 
ing for nutrients during growth, and that the organ nearest to the 
nutrient source will be most successful 37 39 , as nutrients are trans- 
ported to distant organs only after the needs by adjacent organs are 
met 39 . 

As hypothesized, N and P concentrations show different scaling 
relationships among organs. Consistent with Kerkhoff et al. 16 , we 
find that scaling relationships of N concentration in photosynthetic 
organ (leaf) vs. N concentration in non-photosynthetic organs (stem 
and root) tend to be allometric, while that of N concentration 
between stem vs. root is very close to (although not) isometric 
(Table 1). Particularly, there is a faster increase of N concentration 
in non-leaf organs as leaf N concentration increases {b > 1.0), partly 
because non-leaf organs are closer to the nutrient source (soil) 37 39 . A 
second possible reason might be related to the activity of vascular 
tissues. Phloem tissues require more N investments and higher rates 
of N cycling for the high rate of photosynthate transportation and 
photosynthetic activity 16,40 . Among the non- photosynthetic organs, 
we also observe that the scaling relationships are not strictly isomet- 
ric, partly because of the difference in distance to the nutrient 
sources 37-39 . Contrary to N concentration, the scaling slopes of P 
concentration between non-root organs are close to 1, but higher 
in root vs. non-root organs, partly because that root is the closest 
organ to P source. This difference between the scaling relationships 
of N and P concentrations among organs might because of different 
requirements of N and P among plant organs to perform various 
functions. However, it is difficult to fully understand these differ- 
ences without measuring important physiological processes of plant 
organs, such as photosynthesis, respiration, nutrient absorption and 
transportation, or plant functional traits which can represent these 
physiological processes. Future studies including plant physiological 
processes will help to further explain the mechanisms behind the 
scaling relationships of N and P across plant organs. 

In contrast to concentrations, the scaling slopes of N and P 
storages among plant organs are quite similar, which indicates that 
the allocations of N and P in plants largely depends on the allocation 
of biomass. This result is reasonable because both N and P are only 
small portions of total plant tissue, and the pools of N and P in plant 
tissues mainly reflect the tissue size rather than composition. 

Most scaling relationships remain unchanged after phylogeny was 
controlled, indicating that the scaling relationships of nutrient con- 
centrations across plant organs are independent of phylogenetic 
relationships between species. We also observed the same scaling 
relationships for most organ pairs between legume and non-legume 
shrubs. Although plant N and P concentrations are considered phy- 
logenetically conservative traits 16,35,41 , that was not the case for the 
allocation of N and P among organs in this study. The scaling rela- 
tionships of nutrient concentration across most organs are not sim- 
ply a result of shared evolutionary history of these species 16 . 

However, we still found significant difference between species and 
phylogeny level relationships for root N concentration vs. stem N 
concentration. Similarly, we found significantly different scaling re- 
lationships in root N concentration vs. stem N concentration 
between legumes and non-legumes. These results indicate that the 
N fixing pathways used by legumes impact the allocation strategy of 
N in shrubs. The difference in the scaling of N between legumes and 
non-legumes supports the hypothesis that plant functional group 
affects plant traits as well as their scaling relationships 16,42 . There- 
fore, it is important to take plant functional group into consideration 
when exploring plant nutrient allocation strategies. 

We did not observe a significant influence of STN and STP on the 
allocation of nutrients among organs, suggesting that plants do not 
change their nutrients allocation strategies even under soil nutrients 
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Figure 2 | Comparisons of the RMA regression slopes of N concentration and storage (solid dots), P concentration and storage (open cycles) 
among different STN levels (a), (d), STP levels (b), (e) and climatic regions (c), (f). Line bars show 95% confidence intervals (CI). Letters above the line 
bars show the results of likelihood ratio tests (uppercase for N concentration and storage, lowercase for P concentration and storage). Slopes with the same 
letters are not significantly different (P > 0.05), while those with different letters are significantly different (P < 0.05). 



deficiency. Available soil N and P can be absorbed directly by plants, 
and thus might influence N and P concentrations in plants 43,44 . 
However, recent studies on leaf traits of shrubs and grasses con- 
ducted across Inner Mongolia suggests that total and available soil 
nutrient concentrations explain similar and a very small amount of 
variance in leaf N and P concentrations of shrub species 45,46 . These 
results indicate that in the region studied, total soil nutrient concen- 
trations do not differ from available soil nutrient concentrations in 
explaining or predicting plant nutrient concentrations, which is con- 
sistent with our results. Other soil properties such as soil age 47 , soil 
freezing and disturbance 17 may impact nutrient allocation in plants; 
however, these factors are not included in the present study. Further 
studies are needed to examine the effect of soil properties on the 
allocation of nutrient in plants. 

Plants in arid conditions tend to have higher leaf N concentration 
in order to exploit light 25 and increase the water use efficiency at the 
expense of lowering N-use efficiency 28,29 . As expected, increases of N 
concentration in non-photosynthetic organs are even faster in humid 
than in arid regions (Figure 2c). Our result supplements these the- 
ories by considering the relationship of nutrient concentrations 
among plant organs. Although the increase of N concentration in 
non-leaf organs is always faster than the increase in leaf N concen- 
tration in all the three climate regions, leaf N concentration tends to 
be higher at given root or stem N concentration under arid region 
than humid region. The scaling relationships of P concentration 
between non-leaf organs and leaves showed similar variation as N 
concentration among three climate regions. Another possible reason 
might be that water is the limiting factor for vegetation growth in 
shrublands of Northern China and increases in water availability 
promotes plant growth 48 . During growth, more photosynthate and 
nutrients will be transported to non-photosynthetic organs 13 . This 
result indicates that shrub species in Northern China display a strat- 
egy that favors higher nutrient concentration in leaves than in 
non-photosynthetic organs in dry environments. The variation tend- 
encies of RMA slopes of N and P storages are almost the same under 
different aridity conditions, which indicates that aridity can influ- 
ence plant nutrient allocation through variation in plant biomass 



allocation. In semi-arid/semi-humid regions, plants tend to allocate 
more biomass to their root systems than to aboveground organs. 

In summary, we investigated the scaling relationship of N and P in 
leaves, stems and roots based on samples of 126 shrub species from 
172 shrubland communities in Northern China. We found that 
plants have different strategies of allocation for N and P. The scaling 
relationships of N concentration across different plant organs tended 
to be allometric between leaves and non-leaf organs, and isometric 
between non-leaf organs. Whilst the scaling relationships of P con- 
centration tend to be allometric for root vs. stem and for root vs. leaf, 
it tends to be isometric for stem vs. leaf. The scaling relationships of 
both concentrations and storages of N and P among organs are 
affected by aridity, but not by soil nutrient concentration. In arid 
environments, plant tend to have higher nutrient concentration in 
leaves at given root or stem nutrient concentration. Most scaling 
relationships of concentrations and storages of N and P across plant 
organs don't exhibit phylogenetic signal, except that of N concentra- 
tion, which could be attribute to the difference of legumes and non- 
legumes in N utilization. We therefore conclude that consistent 
nutrient scaling relationships among plant organs exist in indepen- 
dently evolving lineages, but factors such as climate and N-fixation 
ability lead to variations in these relationships. 

Methods 

Study site and investigation. This study was carried out at 172 shrubland sites 
extending 46.1 degrees in longitude (86.7-132.8 J E) and 18.7 degrees in latitude (33.7- 
52.4 L, N) in Northern China between July and September (mostly July and August) 
2011 (Figure 3). Three 5 m X 5 m plots were selected at each site to represent the 
natural shrubland communities. We identified all individuals to the species level. In 
each plot, leaf, stem and root biomass were harvested, and all biomass were harvests 
separately for each species. For each shrub species encountered in the community, we 
selected three to five individuals and collected fully expanded leaves, stems and roots 
(mostly coarse roots in the top 30 cm of soil) at each site. Plant samples were oven- 
dried and ground after being transported to the laboratory. In total, we sampled 702 
individuals of 126 shrub species from 71 genera and 33 families. All of the shrub 
species sampled were deciduous except the Pinus tabuleaformis in one study site. 

We collected soil samples using three one-meter-deep pits along the diagonal of 
each plot. For each profile, soil at the depths of 0-10, 10-20, 20-30, 30-50, 50-70 and 
70-100 cm was sampled and soil samples from the same depth were well mixed. Soil 
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Figure 3 | Locations of the sampling sites. The background shows the distribution of shrubland biomes in the northern part of China based on the 
"Vegetation map of the People's Republic of China (1 : 1000000)" 36 . The map was generated using ArcGIS 10.3 (http://www.esri.com/). 



samples were air-dried, had roots removed and were ground to pass through a 100- 
mesh sieve. 

Monthly precipitation and mean monthly temperature was obtained from the 
WorldClim website 49 (available at www.worldclim.org; resolution of 1 km). 

Measurements. The plant and soil samples were analyzed at the Measurement Center 
of the Institute of Botany, Chinese Academy of Sciences. Total nitrogen 
concentrations of soil (STN) and plant samples (leaf, stem and root N concentrations) 
were analyzed using an elemental analyzer (2400 II CHNS; Perkin-Elmer, Boston, 
MA, USA) under 950 C for combustion then reduced to 640°C. Total phosphorus 
concentrations of soil (STP) and plant samples (leaf, stem and root P concentrations) 
were analyzed using the molydate/ascorbic acid method after H 2 S04-H 2 0 2 
digestion 50 . We only used STN and STP at 0-10 cm depth interval for the analyses 
because STN and STP at this depth interval were highly correlated with those at the 
other five depth intervals. 

Scaling of N and P among different organs. A scaling approach, Y = aX b , was used 
to examine the covariation in N and P concentrations because both concentration and 
storage of N and P tend to be log normally distributed 16,51 . The power function can be 
expressed in the form of a linear regression equation after log- transforming, where 
the exponent b is the regression slope and a is the regression intercept. The N and P 
storages of organs were calculated by multiplying the concentrations of N and P of the 
organs by the corresponding biomasses. 

We first applied reduced major axis (RMA) to examine the correlation of nutrients 
among organs. An important advantage of RMA regression compared with ordinary 
least square regression is that RMA minimizes sums of squares in X and Y simulta- 
neously. We assigned the scaling relationship between Y and X as isometric when the 
95% confidence interval (CI) of b contains 1; the scaling relationship is otherwise 
allometric. A b above 1 indicates that, on average, Y increases faster than linearly with 
X; whereas a b below 1 indicates Y increases slower than linearly with X 10 ' 16 . In the 
RMA analyses, we set the organs closer to the top of the plant as the X variable. 

To examine the effects of soil nutrient on the relationship of nutrient concentration 
and storage across organs, we divided all individuals into three equal subgroups based 
on soil nutrient (STN and STP) and compared the exponential slopes among different 
soil and climate levels (STN < 0.9, 0.9 < STN < 2.3 and STN > 2.3 for low, medium 
and high STN, respectively; STP < 0.47, 0.47 < STP < 0.63 and STP > 0.63 for low, 
medium and high STP, respectively). 

We also identified three climatic areas to examine the effects of aridity on the 
scaling relationships of nutrient concentrations and storages across plant organs 
based on the aridity index in growing season (GAI) defined as the ratio of total 
precipitation to potential evapotranspiration from May to October 52 . Potential eva- 
potranspiration was calculated using the Thornthwaite equation 53 . GAI increases 
from arid to humid climate. We divided the range of GAI of the studied region in to 
three groups: arid regions, with GAI < 0.59 (the least one third of the GAI, Arid); 
semi- arid/ semi- humid regions, with GAI between 0.59 and 1.17 (the median one 
third of the GAI, Semi); and humid regions, with GAI > 1.17 (the highest one third of 
the GAI, Humid). 

To test the influence of plant functional groups on the scaling relationship of 
nutrient concentrations and storage across plant organs, we further grouped all 
samples into legumes (samples from the family Fabaceae) and non-legumes (samples 
from other families), as the legume shrubs fix nitrogen whereas the non-legume 
shrubs do not. 



A likelihood ratio test 54 was used to test the heterogeneity between RMA regression 
slopes of different groups, i.e., different soil nutrient level, climatic regions, and 
functional groups. 

Phylogenetic tree and phylogenetic analyses. We also looked for the phylogenetic 
signals of all RMA regression slopes. To do this, we first constructed a phylogenetic 
tree for the 126 species using Phylomatic 55 based on APG III topology 56 . The branch 
lengths were determined using BLADJ algorithm within the Phylocom software 
(http://www.phylodiversity.net/phylocom/) and the node ages were fossil- 
estimated 57 . We then applied a phylo genetically independent contrast (PIC) 
analysis 58 . PIC analysis allows us to calculate N-l standardized contrasts of a trait 
using a series of trait data for N species. We then conducted RMA regression at both 
species-level data and PICs, we tested the effects of evolutionary history on 
covariation of N and P concentrations and storages across organs by comparing the 
RMA regression slopes at species-level and the PIC analysis using the likelihood ratio 
test. 

All analyses were performed using RStudio with the basic, smatr and picante 
packages (http://www.R-project.org/). 
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